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Phase transitions in quark matter and behaviour of physical 
observables in the vicinity of the critical end point 
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Abstract. We study the chiral phase transition at finite T and /is within the framework of the SU(3) 
Nambu-Jona-Lasinio (NJL) model. The QCD critical end point (CEP) and the critical line at finite tem- 
perature and baryonic chemical potential are investigated: the study of physical quantities, such as the 
baryon number susceptibility near the CEP, will provide complementary information concerning the order 
of the phase transition. We also analyze the information provided by the study of the critical exponents 
around the CEP. 
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Understanding the behaviour of matter under extreme 
conditions is one of the most challenging issues in the 
physics of strong interactions. In the last years, major the- 
oretical and experimental efforts have been dedicated to 
the physics of relativistic heavy~ion collisions, looking for 
signatures of the quark gluon plasma (QGP). A region of 
particular interest, that has been so far explored, is the 
T — iiB phase boundary in order to try to explore differ- 
ent regions of the QCD phase diagram. The existence of 
the CEP in QCD was suggested in the end of the eight- 
ies, and its properties have been studied since then (for 
a review see Ref. [T]). The most recent lattice results for 
the study of dynamical QCD with Nf = 2 + 1 staggered 
quarks of physical masses indicate the location of the CEP 
at Te = 162 ± 2MeV, = 360 ± 40MeV [2], however its 
exact location is not yet known (the location of the CEP 
depends strongly of the mass of the strange quark). At 
the CEP the phase transition is of second order, belong- 
ing to the three-dimensional Ising universality class, and 
this kind of phase transitions are characterized by long- 
wavelength fluctuations of the order parameter. 

As pointed out in |3j, the critical region around the 
CEP is not pointlike but has a very rich structure. This 
critical region is defined as the region where the mean field 
theory of phase transitions breaks down and nontrivial 
critical exponents emerge. The size of the critical region 
is important for future searches for the CEP in heavy ion- 
collisions. Some studies have been done in the SU(2) sector 
[3l|4] but less attention has been given to the effects of the 
strange quark. In this paper we aim to investigate the 
chiral phase transition and the CEP in quark matter with 
strange quarks. 



We perform our calculations in the framework of the 
three-flavor NJL model, including the determinantal 't 
Hooft interaction that breaks the {/^(l) symmetry, which 
has the following Lagrangian: 



C — q (id ■ "f — m) q 
+ go det[q{l 



i^Y.[{qX^qf + {q{^J,)X'^qy 
a=0 

75)g] -|-det[^(l -75)q] 



(1) 



Here q ~ (u, d, s) is the quark field with three flavors, 
Nf — 3, and three colors, Nc = 3, m — diag(m„, m^, m-s) 
is the current quark mass matrix and A"^ are the Gell- 

Mann matrices, a = 0, 1, . . . , 8, A" = y^I- The model is 

fixed by the coupling constants gs,9D in the Lagrangian 
([1]), the cutoff parameter A, which regularizes the diver- 
gent integrals, and the current quark masses rui. For our 
numerical calculations we use the parameter set: to„ = 
rud = 5.5 MeV, = 140.7 MeV, gsA^ = 3.67, goA^ = 
-12.36 and A = 602.3 MeV, and we follow the methodol- 
ogy presented in detail in Refs. [5l!6]. 

The baryonic thermodynamic potential has the follow- 
ing form: 



n{n„T) = E-TS- 



(2) 



i—u,d,s 



where i?, S and Ni are, respectively, the internal energy, 
the entropy and the number of particles of the zth quark 
(the expressions are given in Sec. IV of Ref. [5]). The quark 
density for the i quark is determined by the relation 



= NJV = 



+ nt-l)9{A^-p^) (3) 
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Fig. 1. Pliase diagrams in the T — fiB plane. Left panel: go = —12.36. The solid line corresponds to first-order phase boundary 
and the dashed line corresponds to the crossover region. Right panel: go = 0. The lines correspond to first-order phase boundaries 
of the u, d and s sectors (from left). We have used the accepted values ^/ = 10 MeV and ^ly = 12.5 MeV. 



where nf are the Fermi distribution functions of the pos- 
itive (negative) energy state of the ith quark. Along this 
work we impose the condition /ig = so the chemical 
equilibrium condition is /i„ = /i^ = /is = /x^ (with — 

In general, the baryon number susceptibility is the re- 
sponse of the baryon number density psiT^Hi) to an in- 
finitesimal variation of the quark chemical potential fii: 



dfh 



(4) 



This observable is one of the most relevant quantities 
in the context of possible signatures for chiral symmetry 
restoration in the hadron-quark transition and for transi- 
tion from hadronic matter to the QGP [7]. 

The nature of the chiral phase transition in NJL type 
models at finite T and/ or fiB has been discussed by differ- 
ent authors j.5,8j. At zero density and finite temperature 
there is a smooth crossover, at nonzero densities, different 
situations can occur (for details see Ref. [5]). 

The phase diagram for the NJL model is presented in 
Fig. [TJ left panel. For zero temperature the transition is 
of first order. As the temperature increases, the first order 
transition line persists up to a critical endpoint. In the 
CEP the chiral transition becomes of second order. 

We will start our study by investigating the influence of 
the 't Hooft determinant in the T — hb plane in the SU(3) 
NJL model. Once for heavy-ion collisions the isospin chem- 
ical potential (fij) is supposed to be non-zero, the possible 
existence of more than one CEP in the phase diagram was 
proposed in [9]. On the other hand, to study whether a 
kaon condensate is formed, non-zero hypercharge chemi- 
cal potential (fiy) is also important. Having these facts in 
mind, we start by taking gD = and, if we fix nonzero 
values of fii and /iy , there are three first-order phase tran- 
sitions at low T and high /is and thus three CEP as we 
can see in Fig. [1] right panel. However, in our specific 



case, where gu = —12.36, there is only one first order 
transition line and, consequently, only one CEP. This is 
a consequence of the fiavor-mixing effects that cannot be 
neglected in the discussion of the phase diagram [TU] . 

A bound to the size of the critical region around the 
CEP can be found by calculating the baryon number sus- 
ceptibility (xb) and its critical behaviour. If the critical 
region of the CEP is small, it is expected that most of the 
fluctuations associated with the CEP will come from the 
mean field region around the critical region [3]. 

In the left panel of Fig.[2]the baryon number density for 
three different temperatures around the CEP is plotted. 
For temperatures below T'-^^p we have a first order phase 
transition and xb has a discontinuity (right panel of Fig. 

. For T — T'-"^^ the slope of the baryon number density 
tends to infinity at fxs = Ms^^ which implies a diverg- 
ing susceptibility (this behaviour was found in [30 using 
different models in the SU(2) sector). For temperatures 
above T'-^^p ^ in the crossover region, the discontinuity of 
Xb vanishes at the transition line, and the density changes 
gradually in a continuous way as we can see in the right 
panel of Fig. [2l 

We will focus our attention on the critical behaviour 
of Xb in the vicinity of the CEP. At the CEP the baryon 
number susceptibility diverges with a certain critical ex- 
ponent. As pointed out in [TT], the form of the divergence 
depends on the route by which one approaches the crit- 
ical point. If the path chosen is asymptotically parallel 
to the first order transition line, the divergence of the 
baryon number susceptibility scales with an exponent 7^. 
In mean-field approximation, it is expected — \ for this 
path. For any other path not parallel to the first order 
line, the divergence scales with the exponent e = 1 — 1/5. 
Once in the mean-field approximation (5 = 3, we will have 
e = 2/3 and 7, > e is verified. The last condition is re- 
sponsible for the elongation of the critical region, xb being 
enhanced in the direction parallel to the first order tran- 
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Fig. 2. Baryonic density (left panel) and baryon number susceptibility (right panel) as a function of /is for different temperatures 
around the CEP: T^'^^ = 67.7 MeV and T = T^-^^ ± 10 MeV. 
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Fig. 3. The baryon number susceptibility xs as a function of 
Ims - Ms^^l at fixed temperature T'^^^ = 67.7 MeV. 



sition line; we verified that this is also true in the SU(3) 
NJL model [H]. 

For the baryon number susceptibility we will use a path 
parallel to the /i^-axis in the (T, /iB)-plane from lower /i^ 
towards the critical — 318.4 MeV at fixed tempera- 

ture T'^EP = 67.7 MeV. In Fig.[2]we plot xb as a function 
oi \hb ~ f^B^^l close to the CEP. To calculate the critical 
exponent e we will use a linear logarithmic fit 

InxB = -ehi l/iB - ^§^-^1 + const., (5) 

where the term const, is independent of hb- The result 
that we obtain is e = 0.67±0.01, which is consistent with 
the mean field theory prediction: e = 2/3. 

We have analyzed the phase diagram in the SU(3) NJL 
which reproduces the essential features of QCD, such as a 



first order phase transition for low temperatures and the 
existence of the CEP. We verified that, in the chiral limit 
and differently from the SU(2) model, there is no TCP, 
which agrees with what it is expected: the chiral phase 
transition is of second order for Nf = 2 and first order for 
Nf > 3 [13 . A more rich scenario, in SU(3), is found when 
only the nonstrange quark masses are zero (the location 
of the CEP and the TCP depend strongly on the strange 
quark mass) [12]. Around the CEP we have studied the 
baryon number susceptibility which is related with event- 
by-event fiuctuations of /x in heavy-ion collisions. We also 
conclude that in our model the critical exponent e ob- 
tained is consistent with the mean field value e — 2/3. 

Work supported by grant RFBR 06-01-00228 (Yu. Kali- 
novsky) and by grant SFRH/BPD/23252/2005 from FCT 
(P. Costa). 
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